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An introduction and set of guidelines for finite element dynamic modeling of nonrigidized inflatable structures 
is provided. A two-step approach is presented, involving 1) nonlinear static pressurization of the structure and 
updating of the stiffhess matrix and 2) linear normal modes analysis using the updated stiffness. Advantages 
of this approach are that it provides physical realism in modeling of pressure stiffening, and it maintains the 
analytical convenience of a standard linear eigensolution once the stiffness has been modified. Demonstration of 
the approach is accomplished through the creation and test verification of an inflated cylinder model using a large 
commercial finite element code. Good frequency and mode shape comparisons are obtained with test data and 
previous modeling efforts, verifying the accuracy of the technique. Problems encountered in the application of the 
approach, as well as their solutions, are discussed in detail 
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Nomenclature 

element forces in nonlinear static solution, N 

estimates of tangent stiffness, N/mm 

applied loads, N 

unbalanced loads, N 

displacement, mm 

applied load increments, N 

displacement prediction and correction, mm 


Introduction 

I NFLATED cylindrical struts constructed of thin polymer films 
have considerable practical application and potential for use as 
components of inflatable concentrator assemblies, antenna struc- 
tures, and space power systems. Because of their importance, it is 
of great interest to characterize the dynamic behavior of inflatable 
cylinders both experimentally and analytically. It is very helpful to 
take a building-block approach to modeling and understanding in- 
flatable assemblies by first investigating in detail the behavior of 
components such as struts. The film material used for construction 
of such cylinders (and the film-pressurized air interaction) is highly 
nonlinear, with modulus varying as a function of frequency, temper- 
ature, and level of excitation. Thus, the primary motivation of tests 
and analytical modeling efforts was to determine and understand the 
response of inflatable cylinders for different pressures and boundary 
conditions. In this investigation, it is assumed that the structure has 
not yet been (or will not be) rigidized following its inflation. 

In recent years, inflatable structures have been the subject of re- 
newed interest for space applications such as communications an- 
tennas, solar thermal propulsion, and space solar power. 1-6 A major 
advantage of using inflatable structures in space is that they are 
extremely lightweight Additional obvious advantages include on- 
orbit deployability and related volume savings in the launch config- 
uration. A recent technology demonstrator flight for inflatable struc- 
tures was the inflatable antenna experiment that was deployed on 
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orbit from the Shuttle Orbiter. Although difficulty was encountered 
in the inflation/deployment phase, the flight was successful overall 
and provided valuable experience in the use of such structures. 1 * 2 

The Solar Orbital Transfer Vehicle 3 is a potential technology 
demonstrator for solar thermal propulsion. The basic concept be- 
hind solar thermal propulsion is to utilize sunlight or solar energy 
as a means of heating a working fluid (propellant) to provide thrust 
at increased specific impulse. As described in Ref. 7, thrust is pro- 
duced by expanding the heated propellant through a nozzle. No 
combustion occurs, and the thrust level is low. For this reason, solar 
thermal propulsive systems are mainly applicable for orbital transfer 
vehicles. 

Another technology demonstration concept for solar thermal 
propulsion is the Solar Thermal Upper Stage (STUS), which is de- 
scribed in Ref. 8. The engine system envisioned for the STUS is 
designed to utilize hydrogen propellant to produce a thrust level of 
about 8.9 N (2 lbf). Two inflatable parabolic collectors could be used 
that would be rotated and gymballed for focusing sunlight into an 
absorber cavity. The collectors would be inflated after separation of 
the upper stage from the launch vehicle. 

In Figs. 1 and 2, a prototype inflatable solar concentrator 9 * J0 is 
shown that consists of a torus/lens assembly supported by three 
struts. This concentrator is constructed of Kapton^ polyimide film, 
with epoxy as the primary adhesive for joints. Adhesive tape was 
also utilized in this concentrator to attach the lens to the torus. The 
Fresnel lens of such a concentrator assembly would focus sunlight 
into a collector near the fixed ends of the struts. Solar energy stored 
in the collector could be utilized to heat a propellant, as described 
earlier, or for generation of electrical power. The inflatable struts are 
attached to a base plate or spacecraft by means of three cylindrical 
appendages. These hollow appendages also allow inflation of the 
concentrator assembly through air hoses connected at each strut 
It can be seen that inflatable cylindrical struts (Fig. 3) are critical 
components of structural assemblies for practical applications. In 
view of their importance, structural dynamic and static behavior of 
typical inflated film struts need to be investigated. 

The purpose of the current paper is to present a practical finite 
element approach, using a commercially available code, that will 
enable researchers to model the dynamics of nonrigidized inflat- 
able structures accurately. To date, relatively little work has been 
published on dynamic modeling of infiatables, and most methods 
presented are either approximate and lacking physical realism, or 
very laborious in their implementation. First, the paper will give 
an overview of previous modeling of inflatable cylinders. An intro- 
duction is then given for the present method of modeling, which 
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fig. 1 Concept for solar thermal propulsion system utilizing inflatable 
components. 



Fig. 2 Prototype solar concentrator in test chamber. 


involves a two-step procedure for 1) nonlinear pressurization of the 
structure and updating of the stiffness - matrix and 2) modal anal- 
ysis using the updated stiffness matrix. The specific mechanics of 
the method are discussed in parallel with a discussion of their im- 
plementation in NASTRAN code. Problems encountered with the 
method are then discussed, and approaches for their solutions are 
given. The methodology is verified through the correlation between 
an analytical model and modal test results of a thin- film strut 

Recommendations are given for the needed advancement of this 
method and capability to model thin film structures accurately. Fi- 
nally, conclusions are drawn on the usefulness of the methodology. 

Review of Previous Work with Inflated 
Beamlike Structures 

A number of investigators have considered the use of inflatable 
cylindrical beams. Perhaps the earliest was Otto, 11 who published 
ideas for inflated tubular frames for use in structures such as orbiting 
platforms. A more recent proposed application involves the use of^ 
inflatable beam segments to replace solid segments of the^ace^ 
^buttle remote manipulator system and, thus, reduce storage space ^ 
and inertia of the arm. 12 

Several papers on static structural analysis of inflated cylinders 
describe different techniques such as linear shell theory and non- 
linear and variational methods. 1 3 “ 17 Fichter 18 developed nonlin- 



Fig. 3 Inflatable cylindrical strut in cantilever test configuration. 


with internal pressure and axial force set to zero, reduced to the 
Timoshenko beam equations. Comer and Levy, 19 Webber, 20 and 
Main et al. 21 modeled an inflated beam using an Euler-Bemoulli 
approach, but with modified material longitudinal modulus and 
beam moment of inertia. Wrinkling behavior of the material was ac- 
counted for by the assumption of no load-carrying capability when 
longitudinal stress in the beam cross section reached zero. 

Additional work of significance involves rigidization of inflated 
beam structures. One proposed concept is the use of injected foam 
that fills the cylindrical beam cross section, subsequently hard- 
ens, and, thus, rigidizes the structure. This approach is discussed 
in Ref. 22. 
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could be utilized in approximating lower-order frequencies of in- 
flatable beams. Main et al. wrote a very significant paper describing 
results of modal tests of inflated cantilever beams and the determina- 
tion of effective material properties. 24 Changes in material proper- 
ties for different pressures were also discussed, and the beam model 
was used in a more complex structure. The paper demonstrated that 
conventional finite element analysis packages could be very useful 
in the analysis of complex inflatable structures. 

Finally, Refs. 25-27 describe closed-form and finite element 
beam representations of inflatable cylinder dynamics, along with 
shell-element models for comparison. Jt was found, that frequency-, 
dependent modulus, or bending stiffness, was required for the beam- 
element models to represent the dynamic test results adequately. 
Shell models performed reasonably well with linear material prop- 
erties. Note, however, that the film modulus was varied in modal 
analyses of the shell-element models of Refs. 25-27 to approxi- 
mate the effects of pressure and still allow for linear eigensolutions. 
Pressure was included directly only in frequency response analy- 
ses. A more realistic and satisfying approach is to include pressure 
loading directly in the modal analysis process. Such an approach is 
investigated in this paper and is described in the following sections. 

Methodology Overview 

The purpose of this investigation was to verify a procedure for 
dynamic modal analysis of nonrigidized inflatable structures in 
NASTRAN that accounts for the increase of structural stiffness due 
to internal pressure. In this study, the MacNeal-Schwendler Cor : 
poration (MSC) version of the finite element code was utilized. 2 * 
Ordinary eigenvalue analysis in MSG/NASTRAN does not allow 
for the effects of loads applied on a model because modal char- 
acterization of a structure is by definition free from active loads. 
However, in the case of inflated (nonrigidized) structures, bending 
stiffness is mainly due to internal pressure. Therefore, a standard 
modal analysis that does not account for internal pressure loads will 
give incorrect results. To solve this problem, the proper stiffness 
matrix for the pressurized structure can be generated in a nonlinear 
static analysis and then imported into a linear, eigenvalue analysis. 
Of course, modeling of rigidized inflatables is considerably more 
straightforward and may not require a nonlinear approach. 

NASTRAN Nonlinear Static Analysis 

The key to this methodology, and the most difficult step, is the 
nonlinear static analysis. To ^understand more fully the problems 
that arise during this procedure, and the reasons their solutions are' 
successful, an understanding of the nonlinear static analysis proce- 
dure (SOL 1 06) in NASTRAN is necessary. This procedure provides 
the ability to handle geometric nonlinearities, kinematic nonlinear- 
ides, follower forces, large displacements, displacement-variable 
loads, material nonlinearity, nonlinear stiffness relationships, buck- 
ling analysis, contact analysis, and variable boundary conditions. In 
this paper, only the characteristics of the procedure that are appli- 
cable to modeling of inflatable structures will be discussed. These 
capabilities are the stiffness updating and follower forces schemes, 
including their iteration and convergence procedures. 

The nonlinear static analysis utilized in this study involved the 
updating of the stiffness matrix throughout the procedure over a 
nonlinear curve. NASTRAN is able to do this through its load- 
displacement iteration process, shown conceptually in Fig. 4. This 
process involves initially applying a predetermined fraction of the 
total load on the structure. Through application of this load, and 
the initial geometry and property values of the loaded elements, ini- 
tial values for the displacement and elemental forces are obtained. 
When the element nodal (edge) forces and the applied load are 
not in equilibrium, the iteration is repeated with updated values for 
the displacement and element forces. These values are calculated 
through the change in displacement due to the difference between 
the applied load and elemental load. This cycle is repeated until 
equilibrium is achieved. The point at which equilibrium occurs is 
called the convergence point It is at this point that the stiffness 
matrix is updated from the last iterated displacement and applied 
load. After convergence of the initial applied load, the analysis ad- 
vances and applies a greater fraction of the combined load. A solu- 
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Fig. 4 Basic concept of nonlinear static solution procedure utilized in 
MSC/NASTRAN. (Copyright © 1994 by MacNeal-Schwendler Corpo- 
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by permission.) 



tion is obtained when convergence occurs for the entire load. The 
NASTRAN code has several stiffness updating schemes, iteration 
and advancing schemes, and convergence criteria, 2 * and the process 
illustrated in Fig. 4 can vary considerably depending on the choice 
of solution parameters. The manipulation and use of these analysis 
schemes and parameters present the greatest challenge in use of this 
methodology. 

For the thin-film cylinders investigated in this study, convergence 
of the NASTRAN nonlinear pressurization solution was difficult to 
obtain because the thin material tended to have large displacements 
in the initial iterations of the analysis. It was found that if the model 
could be stabilized through the first several iterations (to reach the 
point where changes in nodal displacements between analysis steps 
became smaller) the solution proceeded easily to convergence. The 
authors found that one way to stabilize the nonlinear solution was to 


constrain the models artificially. For example, artificial constraints 
were placed at points on the film surface along the length of a free- 


free pressurized cylinder (Fig. 5), at a sufficient number of locations^ 
to suppress the largest displacements. Only the rotational degrees 


of freedom were allowed at these constrained nodes. This method 


was found to be workable for the analysis, an^yn generator this J?* 
methodology, particularly for large structures, points on the thin- 
film surface must either be artificially constrained as illustrated in 
Fig. 5, or artificially stiffened in the first analysis steps. Models for 
large complex structures may require a trial- and-error approach of 
observing the locations of largest displacements in the initial steps 
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An alternative to the use of artificial constraints is the applica- 
tion of artificial stiffness to the model during the initial load itera- 
tions. This can be done in automated fashion in MSC/NASTRAN 
through the KDIAG parameter, 28 which results in the addition of 
a user-prescribed stiffness value (value of KDIAG) to the diago- 
nal terms of the model stiffness matrix. The additional stiffness is 
applied only during the earliest iterations of the procedure to get 
the nonlinear solver started and is then removed automatically by 
NASTRAN once the model is stabilized. Obviously, the automated 
nature of this approach makes it attractive, and it is the method 
— ^recommended by MSC.-However, there appears to be a potentially 
significant disadvantage to the KDIAG approach for very flexible 
thin-walled structures. The procedure was not utilized in this study, 
but it was discovered during modeling of a large inflatable optical 
membrane that the analysis results were sensitive to the value of 
KDIAG chosen by the analyst As a result, trial and error assess- 
ments of the value to be added to the diagonal stiffness terms may 
be necessary, thus, largely negating the advantages of the KDIAG 
approach. 


Fig. 6 Example deformed shape of inflated strut with artificial con- 
straints. 

of the pressurization analysis and then applying artificial constraints 
at those locations. 

It was determined for the cylinder models that the artificial con- 
straints only marginally changed the final values of the updated 
stiffness matrix, even though the constraints modified the displaced 
shape of the example pressurized model discussed in the preced- 
ing paragraph (Fig. 6). The spikes evident on the surface of the 
model result from the tendency of the film material to have large 
displacements during pressurization, combined with the effects of 
the artificial constraints. Apparently, the stiffness of the thin-film el- 
ements is controlled largely by the magnitude of the displacements 
and not by the direction of them. Even though the artificial con- 
straints altered the geometry of the elements, the strain and stress of 
those elements remained mostly unchanged and the stiffness matrix 
was only slightly altered. 

The minimal effect of artificial surface constraints on model stiff- 
ness can be further explained by noting that the magnitudes of the 
elemental normal forces rely greatly on the curvature of that ele- 
ment. The greater the angle of curvature, the larger the component 
of the tensile force that can react against the applied load. Left 
alone, an element will displace rather than increase its curvature, 
causing convergence problems in NASTRAN due to large displace- 
ments. The artificial constraints allowed the elements to increase 
their rate of curvature in relation to the load applied, rather than to 
displace as a rigid body. This increase in curvature is still translated 
into strain, though in a different direction, allowing stress and, thus, 
updated stiffness to remain mostly unaltered in comparison to an 
unconstrained pressurization. Of course, if two adjacent nodes on 
an element were artificially constrained, the stiffness matrix would 
be significantly modified. 

Once the artificial constraints are in place, the nonlinear static 
analysis is performed. In this analysis, NASTRAN solutioi^OL*106 
allows for the manipulation of the parameters and solutioivschemes 
in the model description (bulk data section). 28 Many convergence 
problems can be overcome through the manipulation of these pa- 
rameters. The authors found that for this study it was adequate to 
alter only the convergence criteria and their margins of error and to 
use default values for other analysis parameters. Further investiga- 
tion into the manipulation of other parameters is certainly needed. 
NASTRAN has the option of choosing between displacement, en- 
ergy, or load as the criterion from which it will base its convergence. 
The authors found it to be useful to activate all three to allow anal- 
ysis advancement from whichever criterion converges first It was 
also found useful to increase the error tolerances for these criteria 
from their default values, though care must be taken in doing so. 
Solution accuracy is, thus, traded off with the increased potential for 
convergence. It was also discovered that too large an error tolerance 
was just as fatal as one too small for achieving conveigence. It is 
thought that too large a tolerance takes later iterations away from 
the conveigence point. 


NASTRAN Modal Restart 

On convergence of the nonlinear static analysis, a modal restart 
is performed (Fig. 7). This restart is simply the second stage of the 
analysis procedure, which is activated through the RESTART com- 
mand in the executive control section 28 of the file. This command 
imports the MASTER file (containing information about the model 
and results of the first run) from the nonlinear analysis and claims 
it as its own. From this MASTER file, the restart uses the original 
geometry of the model and the stiffness matrix from a chosen itera- 
tion. In the bulk data section of the restart run, the desired iteration 
of the nonlinear static analysis stiffness update can be chosen. Also, 
at this point in the analysis, the artificial constraints discussed ear- 
lier must be removed by deleting the appropriate lines in the sorted 
bulk data from the nonlinear analysis. The procedure for removing 
the constraints is demonstrated in the last several lines of Fig. 7, 
where the appropriate line numbers to be deleted were found in the 
sorted bulk data echo in the nonlinear f06 (results) file. Remaining 


$ Restart Program for modes analysis using the 
$ stiffness matrix of a static analysis MASTER 
ASSIGN 0UTPIIT2 - * case23.op2 ‘ , UNIT - 12 
$ Direct Text Input for File Management Section 
$ Normal Nodes Analysis, Database 
$ 

RESTART UERSION-1 ,KEEP . 

ASSIGN MASTER-' case22. MASTER* 



SOL 103 
TIME 6 GO 

$ Direct Text Input for Executive Control 
CEND 

SEALL - ALL 
SUPER - ALL 

TITLE - MSC/NASTRAN job MODES RESTART 
ECHO - NONE 
HAXLINES - 999999999 

$ Direct Text Input for Global Case Control Data 
SUBCASE 1| 

$ Subcase name : rigid 
SUBTITLE-rigid 
METHOD - 10 
SPC - 2 

UECT0R(S0RT1 V REAL)-ALL 
SPCF0RCES(S0RT1 ,REAL)-ALL 
BEGIN BULK 
PARAH,NML00P,12 
EIGRL.10,0.,500.,20 
/,2350 

7 .2351 

7.2352 
/ v 2353 
7,2354 
/,2355 
$ 

ENDDATA 


Fig. 7 NASTRAN runstrcam for modal solution utilizing updated 
stiffness matrix from nonlinear pressure solution. 
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sections of the restart are organized as an ordinary SOL 103 eigen- 
value analysis. Within the restart command, it is useful to use the 
KEEP command for retaining the old MASTER file for later use. 
Reference 28 provides an overview of the capability and file struc- 
ture for the NASTRAN finite element code, and additional refer- 
ences provide further information for dynamic modeling 29 and non- 
linear analysis. 30 For the interested reader, the various commands 
and entries for the restart file in Fig. 7 are described in detail in 
Ref. 28. 


Application of Methodology to Dynamics 
of Inflated Cylindrical Beams 

The solution procedure described in the preceding two sections 
was performed and ultimately verified through the modeling of thin- 
film inflatable struts characterized earlier. A typical inflated cylin- 
drical beam or strut of the type investigated in this paper is shown in 
Fig. 3. The beams studied are 2.44 m (8 ft) long with 1 5.2-cm (6-in.) 
diameter, are constructed of Kapton polyimide film, and are sealed 
at the ends using styrofoam plugs. Openings in the plugs provided 
the means for inserting air hoses to inflate the structure to desired 
pressure. Epoxy was used to bond the foam plug to the polyimide 
film walls of the cylinder. 

The beams were constructed by overlapping the edges of the 
polyimide sheet and placing a thin layer of epoxy adhesive along 
the overlapped area to form a joint running the length of the beam. In 
a cross-sectional view, the bonded joint looks like a sandwich, with 
the epoxy layer between two layers of Kapton film. Of course, the 
bonded region has much higher stiffness than a nominal polyimide 
section of cylinder wall, and this had to be accounted for in modeling 
of the struts. The stiffness of the joint appears to be dominated by 
the epoxy layer rather than the polyimide material, partially due to 
the greater thickness of epoxy. Polyimide film of 51 fxm (0.002 in. 
or 2 mil) thickness was utilized, and the epoxy layer was typically 
255 iim (0.01 in. or 10 mjl) thick. 

Quadrilateral shell elements were used to model the strut film 
material. Nonstructural masses were added to these elements to 
compensate for the mass of air within the structure. A single line 
of elements along the length of the strut was thickened to 293 fim 
(11.5 mil) to compensate for the overlap and adhesive along the 
seam. Solid elements were used for the end plugs. 

Once the strut model was completed, the first step in the modal 
analysis procedure was to constrain the model (artificially in the case 
of free-free boundary conditions) and perform a nonlinear static 
pressurization using NASTRAN SOL 106. The constraints were 
placed on the relatively stiff foam end plug because both cantilever 
and free-free configurations were to be investigated. Results of the 
nonlinear analysis appear very reasonable, as shown in Fig. 8. The 
effect of the seam, with its higher stiffness due to the film overlap 
and adhesive, is observed in the “hotdog” shape of the strut with 
the free end bending upward. During the SOL 106 run, the model 
stiffness matrix was updated to capture the pressure stiffening effect 


The second step in the procedure was to run a linear modal anal- 
ysis using the updated stiffness matrix. For the free-free cases, the 
artificial constraints on the end plug were removed before initiat- 
ing the eigensolution. Figure 7 shows the NASTRAN modal restart 
runstream used for the free-free configurations of the cylindrical 
beam. It is noted that, for this analysis, artificial constraints were 
not placed on the film surface, as was done for the example in Figs. 5 
and 6. However, for large inflatable structures, it appears to be nec- 
essary to constrain the film at points of largest displacement during 
pressurization. 

. The results of analysis using the present methodology compared 
to modal test data 31 and previous modeling results from Refs. 25-27 
can be seen in Table 17 where frequencies for the lowest bending 
modes are shown for each case. For all frequencies in Table 1, the 
mode shapes were classical beam bending- type modes. For exam- 
ple, Fig. 9 shows the second mode shape for the free-free conditions 
and 9.7-kPa (1.4-psi) pressure. Free-free and cantilever configura- 
tions were tested and modeled at multiple pressure levels for each 
boundary condition. As can be seen in Table 1 , the free-free results 
fell within an acceptable margin of 6% frequency error in compar- 
ison to test frequencies. Furthermore, the free-free analytical fre- 
quencies obtained using the present method are considerably more, 
accurate than the results from Refs. 25-27. Results of two models 
from Refs. 25-27 are given in Table 1 , a linear shell-element model 
in which the effects of pressure were represented by increasing the 
film modulus and a nonlinear shell-element model with frequency- 
dependent modulus and directly applied pressure. For the nonlin- 
ear shell model, the modal data were extracted from NASTRAN 
frequency response results in a rather laborious process. 
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Fig. 8 Deformed shape of pressurized strut with end constraints. 


Table 1 Comparison of analytical and experimental modal results for free-free and cantilever boundary 

conditions and various internal pressures 



✓ 

Present method 

Refs. 25-27 linear 

Refs. 25-27 nonlinear 

Mode 

Experiment frequency, Hz 

Frequency, Hz % error Frequency, Hz 

% error 

Frequency, Hz 

% error 

, 

24.30 

22.85 

Free-Free, 5.2 kPa 
5.97 

22.93 

5.64 

26.50 

-9.05 

2 

71.80 

67.83 

5.53 

80.78 

-12.51 

67.40 

6.13 

, 

24.60 

25.30 

Free-Free, 9.7 kPa 
-2.85 

23.05 

6.30 

26.20 

-6.50 

2 

71.40 

70.10 

1.82 

81.20 

-13.73 

68.90 

3.50 

1 

2.78 

3.61 

Cantilever, 3.5 kPa 
-29.86 

2.18 

21.58 

2.90 

-4.32 

2 

25.33 

26.69 

-5.37 

28.78 

-13.62 

24.10 

4.86 

3 

78.73 

75.00 

4.74 

85.83 

-9.02 

67.20 

14.64' 

1 

2.81 

4.44 

Cantilever, 6.9 kPa 
-58.01 

2.19 

22.06 

2.90 

-3.20 

2 

25.10 

28.22 

-12.43 

28.91 

-15.18 

23.90 

4.78 

3 

79.04 

76.58 

3.11 

86.21 

-9.07 

67.10 

15.11 
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Fig. 9 Second bending mode for free-free inflated strut, pres- 
sure = 9.7 kPa (1.4 psi). 


Results for the present methodology in Table 1 for the cantilever 
arrangement show good model-to-test agreement for the second and 
third bending modes, but large frequency errors are seen for the fun- 
damental modes. It is possible that these errors in the fundamental 
frequencies are due to film wrinkling effects near the fixed end of 
the strut that were not taken into account in the finite element model. 
Important studies of membrane wrinkling behavior, including can- 
tilever inflated cylinder examples, are described in Refs. 32 and 33. 
Errors in higher-order mode shapes are possibly due to nonlinear 
elastic modulus as a function of frequency. Note that geometric non- 
linearity (film stretching and stiffening due to internal pressure) was 
taken into account in modeling efforts described in this paper, but 
the use of nonlinear material properties in the present methodology 
also needs to be investigated. With the exception of the fundamental 
cantilever frequencies, better overall accuracy was achieved with 
the nonlinear pressurization and linear modal approach than was 
obtained in Refs. 25-27. 

Numerical accuracy of the finite element modeling procedure in- 
vestigated in this study was verified not only by comparison to test 
results, but also by comparison to previous finite element modeling 
efforts and Euler-Bemoulli beam solutions (for the same inflat- 
able strut described in this paper) documented in Refs. 25-27. Fur- 
thermore, and perhaps more importantly, because cylindrical beams 
were modeled, the bending modes were expected and observed to 
be classical free-free and cantilever beam modes. Prediction of the 
correct mode shapes by the models developed for this study, as illus- 
trated by Fig. 9, provided tremendous overall confidence in the nu- 
merical procedure. Thus, it becomes a matter of verifying numerical 
accuracy for the model frequencies; Shell models in Ref. 26 had ap- 
proximately 2000 quadrilateral elements. In the current study, fewer 
quadrilateral elements (832) were utilized, but it can be concluded 
from the comparison with Refs. 25-27 and test data in Table 1 , 
particularly for frequencies of the free-free cases, that a sufficient 
number of elements were used. 

Where the use of artificial constraints in the NASTRAN non- 
linear static solution is concerned, it was found in this study and 
related investigations 34 * 35 that the locations chosen and the degrees 
of freedom constrained could be varied without significant effect on 
the frequencies and mode shapes. Stress distributions and overall 
displaced shapes obtained in the nonlinear analysis also appeared 
generally reasonable. 

Accuracy of test data was ensured through a number of 
procedures, 31 including calibration assessments for the laser vi- 
brometer, load cell, impact hammer, and signal conditioner; func- 
tional testing of the test setup; linearity checks for response func- 
tions; and visual inspection of instrumentation for verifying proper 
locations and orientation. Comparison of test frequencies for dif- 
ferent pressure levels in the struts also provides a valuable means 
for assuring consistency and accuracy in the modal data. Test im- 
plementation of the desired cantilever boundary conditions was ac- 
complished through the use of a metal clamp attached to the foam 


plug at the upper end of the strut (Fig. 3). Free-free conditions were 
achieved through the use of a flexible rubber cord at the upper end. 
More detailed information on the test configurations and procedures 
can be found in Refs. 26 and 31. 

Note that the analytical modal results in Table 1 were obtained in 
the finite element models without the manipulation of the original 
element properties. That is, the manufacturer's modulus, density, 
and Poisson’s ratio data were used for the polyimide film, which 
provides a high level of confidence in the solution approach and 
its ability to predict dynamic characteristics of inflatable structures. 
- This methodology has also been implemented successfully in more 
sophisticated NASTRAN models, which are described in Refs. 34 
and 35. Positive results from.these larger and more complex models 
furthered the confidence of the authors in the usefulness and ver- 
satility of the procedure. Model-to-test comparisons shown in this 
paper are only for bending motions of the struts and do not involve 
the shell behavior. Analytical and experimental comparison of shell 
modes was not made because the modal test procedure was limited 
in its ability to identify shell modes. 

Important and independent work involving dynamic modeling 
of an inflatable torus was done in approximately the same time 
frame as this effort. 36 The same general approach presented in this 
paper (static pressurization and subsequent modal analysis using the 
updated stiffness) was used in Ref. 36 with good results, providing 
further verification that this approach is proper and accurate for 
modeling thin-film inflatable structures. 

Advantages and Limitations of the Procedure 
for Modeling Thin-Film Inflatable Structures 

An obvious major advantage of the approach described in this 
paper is the physical realism achieved by accounting for the internal 
pressure directly in modal analysis and by accurately modeling the 
film stiffening due to pressure. Also, this methodology provides the 
convenience of linear modal analysis in NASTRAN as the second 
step in the process. The work described in this paper is an extension 
of previous research described in Refs. 25-27, where for NASTRAN 
modal solutions the pressure stiffening effect was included indirectly 
by modifying the strut film modulus. Pressure was modeled directly 
in Refs. 25-27 only for frequency response analysis incorporating 
nonlinear modulus, and modal properties were then obtained in a 
rather painstaking effort from the response functions. The current 
approach is considerably more convenient for the analyst. 

Apparent limitations in this methodology were found throughout 
the course of the investigation. The small thickness of the inflatable 
structure material is apparently the primary factor causing difficulty 
in the use of this methodology. In simple terms, very thin inflat- 
able structures rely on internal pressure for bending stiffness of the 
membrane. The smaller the thickness of the pressurized structure, 
the larger the nodal displacements in the initial load steps, making 
it difficult to obtain a stable analysis. It is also known that elemen- 
tal forces calculated and applied during the NASTRAN iteration 
process are a function of the volume of each element Decreasing 
the thickness of the film, and, thus, the volume of the elements 
and the elemental forces applied in each load step, causes a greater 
number of iterations to be needed for convergence, decreasing the 
likelihood of NASTRAN reaching a solution. Parameter studies for 
large inflatable structures have shown that increasing the material 
thickness generally eliminates the convergence problems discussed 
in this section. 

The angles between very thin adjacent elements also affect the 
convergence of the nonlinear pressurization solution. This is an im- 
portant factor, because the element tensile forces, which balance 
the pressure loads and are proportional to the angles between ad- 
jacent elements, impart transverse stiffness to the membrane and, 
thus, stabilize the numerical solution. The initial element geome- 
tries, of course, do not have curvature. Thus, the initial elemental 
forces that balance pressure loads are not due to the curvature, but 
rather the angle between adjacent elements. Examining a group of 
elements as a free-body diagram (Fig. 10), it becomes apparent that 
the forces reacting to the applied pressure load come from the com- 
ponents of the element edge forces that are parallel and opposite to 
the pressure vector. These forces are applied to the center element 
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Fig. 10 Quadrilateral plate element deformation under pressure 
loading. 


nodes through the surrounding elements. The smaller the angle be- 
tween adjacent elements, the smaller the force components to resist 
the pressure, and, thus, initial displacements become larger and re- 
tard convergence of the solution. Two difficulties related to small 
angles between adjacent elements are discussed in the following 
paragraphs. 

The use of a very fine mesh was found to affect convergence ad- 
versely in this particular analysis. Not only does a finer mesh com- 
plicate the model and increase the computational run time, it also 
decreases the angles between elements along the curvature of the 
film. Therefore, a greater number of artificial constraints are needed 
to provide equilibrium for the elements in early stages of loading. 
Because of this phenomenon, a coarser mesh may be advantageous, 
though this observation has not been proven in general for very thin 
pressurized structures. The fineness of the mesh may not be an is- 
sue if the artificial applied stiffness approach (KDIAG parameter 
discussed in the “Methodology Overview” section) is utilized. 

Greater difficulty has been encountered in obtaining numerically 
stable solutions for geometrically larger thin-film structures 35 com- 
pared to smaller structures. It appears that important parameters are 
the ratios of material thickness to major dimensions such as length 
or diameter. Decreases in the angles between adjacent elements and 
an increase in the/nuraber of elements needed to characterize the 
structure properly are results of a geometrically larger model. As 
described in this section, both of these conditions can decrease the 
ability of NASTRAN to converge on a solution in its nonlinear 
analysis. 

Finally, note that NASTRAN is a large, general-purpose finite 
element code that is not optimized for nonlinear analysis. The re- 
quired use of artificial constraints or artificial stiffness as described 
in this paper clearly points out the difficulty that can be encountered 
in using general-purpose finite element codes for inflatable struc- 
tures modeling. Finite element packages more suitable for nonlin- 
ear analysis should be investigated, or, if necessary, developed, to 
make the modeling of thin inflatable structures more convenient for 
analysts. However, it has been clearly shown in this study that accu- 
rate dynamic results can be obtained for inflatable structures using 
commercially available finite element codes. 

Summary 

The investigation and verification of a two-stage finite element 
procedure for determining the natural mode$ and frequencies of very 
thin inflatable structures has been described in this paper. Geomet- 
ric nonlinear static analysis, which accounts directly for pressure 
stiffening of the model and results in an updated stiffness matrix, is 
the first step in the procedure. Subsequently, a linear modal analysis 
is conducted using the modified stiffness matrix from the first step. 
Obvious advantages of this approach are that it provides physical 
realism in modeling of pressure stiffening and that it maintains the 
analytical convenience of a standard eigensolution once the stiffness 
has been modified. 


Though not investigated in this paper, materia] nonlinearity can 
also be accommodated in the first stage of this approach, and the 
subsequent linear modal analysis can be done, avoiding a painstak- 
ing process of determining the modal properties from a nonlin- 
ear NASTRAN frequency response analysis. Although good results 
have been shown in this paper using only the geometric nonlin- 
ear capability of the finite element code, use of nonlinear material 
properties in this procedure should be a subject of future research. 

The purpose of this investigation was to verify the usefulness of 
the methodology through the creation of a working and accurate 
dynamic-model-of-an-inflated thin film cylinder, and this has been 
accomplished. This work has contributed to the understanding of 
appropriate finite element modeling procedures for inflatable struc- 
tures, and should help to advance the technology of inflatable space 
systems. 
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